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Abstract 

We present a model that describes states of photon pairs, which have been generated by biexciton 
cascade decays of self- assembled quantum dots, use of which yields finding that agrees well with 
experimental result. Furthermore, we calculate the concurrence and determine the temperature 
behavior associated with so-called entanglement sudden death, that prevents quantum dots emit- 
ting entangled photon pairs at raised temperatures. The relationship between the fine structure 
splitting and the sudden death temperature is also provided. 
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Entanglement, the intriguing correlations of quantum systems, has found extensive appli- 
cation in the areas of quantum computation [l| and quantum communication j^j. Although 
this fundamental feature has attracted a lot of attention since the early years of quantum 
mechanics, it is still far from being completely understood. In recent times, experimenta- 
tion has found that entanglement within a two-qubit system can disappear abruptly within 
a finite time during the decoherence process [3H5|. This phenomenon, termed "entangle- 
ment sudden death" (ESDJ^ has been demonstrated in optical systems under rather special 
artificial environments 

The generation of entangled photon pairs by using a biexciton cascade process in a single 
quantum dot (QD) was first proposed by Benson et al. |8], and has been recently realized 



by several groups [9, In this process, two electrons and two holes are initially generated 
in the QD (biexciton XX), and then a biexciton photon Hxx or Vxx is emitted as the dot 
decays to an exciton (X) state by recombining one electron and one hole. The polarization 
of the biexciton photon is either horizontal (H) or vertical (V"), in accord with the decay into 
the exciton state Xh or Xy, respectively. After a time delay r, the other electron and hole 
recombine to emit an exciton photon H x or Vx with the same polarization as that of the 
earlier biexciton photon. When the two recombination paths (H or V^)are indistinguishable 
in frequency, the output is the maximum entangled Bell state |$) = (\HH) + \VV))/y/2. 

Unfortunately, in QDs, due to the underlying atomic symmetry being C% v rather than 
Ca v Uj, Il2|], the intermediate two exciton state has a small splitting of about several tens 



of /xeV, much larger than the intrinsic width of the emission line (~ //eV). This splitting, 
labeled as S in Figure 1, has conventionally been called the "fine structure splitting" (FSS). 
Due to this FSS, the actual output is not a pure polarization entangled state, but a state 
with both polarization and frequency entangled, along with a phase difference between H 
and V Q, 

|*> = -j=(\H xx H x ) + e iST \V xx V x )). (1) 

If the frequency is traced over, this state becomes just a mixed state with classical correla- 
tions; only when S ~ can we achieve an entangled state. 

In experiments, there are many methods to tune the 
ing uniaxial stress 15] , or applying a magnetic field [16 
various other experiments are designed to erase the frequency information by filtering the 
photons spectrally [s| or selecting the exciton (X) photon within only a small emission delay 



'SS, such as by annealing [yj], exert- 
or an electric field [n]]. Meanwhile, 



relative to the biexciton (XX) photon [13], by which entanglement is greatly enhanced. 

In this work, taking the acoustic phonon-assisted transition process into account, we 
investigated entanglement for photon pairs generated from a QD. A similar ESD, which 
depends on temperature and FSS, has been found in this process. 

A schematic of the biexciton cascade process in a QD is presented in Fig.[TJ where the 
vacuum state (G), the two intermediate exciton states (Xy and Xh), and the biexciton state 
(XX) are labeled respectively as |0) ~ |3) for convenience. 

The QD is initially excited to the biexciton state by a short-pulsed laser, and subsequently 
evolves freely. Considering the phonon assisted process, the density matrix (p) of this four- 
level system can be derived using the master equation (h = 1) jis| . 



p = -i[H ,p] +L(p), 



(2) 



where 



H = 

i=0 

The dissipation term in the Lindblad form is described by 



(3) 



L(P) = -[ 7 3 2 (|2)(3|p|3)(2|-|3)(3|p-p|3)(3|) 
+73i(|l)(3|p|3)(l|-|3)(3|p-p|3)(3|) 
+72o(|0)(2|p|2)(0|-|2)(2|p-p|2)(2|) 
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FIG. 1: Energy level schematic of the biexciton cascade process, the ground state G (|0)), the 
two linear polarized exciton state Xh (|2)) and Xy (|1)) and the biexciton state XX (|3)). The 
spontaneous emission process is denoted by 7^ (712 and 721 are phonon assisted transition rates). 
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+7io(|0)(l|p|l)(0|-|l)(l|p-p|l)(l|) 
+ 721 (|l)(2|p|2)(l|-|2)(2|p-p|2)(2|) 
+ 7l2 (|2)(l|p|l)(2|-|l)(l|p-p|l)(l|)] 



(4) 



where 732, 731, 720, 7io are the spontaneous emission rates, while 721 and 712 denote the 
phonon assisted transition rates between |1) and |2). The phonon absorption rate of the 
state |1) is 712 = kNb and the emission rate of |2) is 721 = k(Nb + 1), where n is the 
phonon-QD interaction rate, which is approximately proportional to the cube of the energy 
splitting S [19], and N B represents the Bose distribution function of a phonon with energy 
S, N B = [exp(S/k B T) - I}' 1 . 

This emission process produces radiation in a mixed state, from which the two-photon 
coincidence measurements single out the polarization density matrix of the photon pair 
Ppoi, which is projected onto the subspace spanned by the four basis {\HiH 2 ), \H1V2), 
\V1H2), IV1V2)}. Its matrix elements, experimentally reconstructed by means of quantum 
tomography, theoretically read as [20 ] 



(Pi^l Ppoi IGC2) = A dt dt' 



tm t r 
rt 



x (<(*R(*>c^':m*)) 



(5) 



where (p, u, £, £) {H,V}, and the as are the dipole transition operators, a H = |3)(2|, 
<j v = |3)(1|, o = |2)(0|, a V2 = |1)(0|. A is a normalization factor to ensure that the trace 
of pp i is unitary The limits t m and tu (t' m and t' M ) defined the temporal window related 
to the detection of the photons from the biexciton cascade process. The emission delay of 



the exciton photon relative to the biexciton photon, see in Eq. (JT]), is r = t' — t. A 



correlation functions are computed by applying the quantum regression theorem 21]. 
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Ideally, the polarization density matrix calculated from Eqs. (j2J) - (jSJ) has the form 

( pn p u \ 
p 22 

p 33 
V p 4 l p 44 ) 

All non-diagonal terms except p 14 and P41 are zero, and the physical reason for this is that 



Ppol 



(6) 
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the Hamiltonian in Eq.(j2]) only couples states that share the same excitation labels H and 
V. 

In realistic situations, aside from phonon-induced spin scattering, the system may be 
affected by some other decoherence processes, such as decay-path distinguishability, back- 
ground noise, cross dephasing and pure dephasing. However, cross- dephased light has been 
shown to be weak [22], and photon pair emission in a QD is robust against single photon 
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221 ] . Thus we can ignore cross 



decoherence being typically limited by pure dephasing 
and pure dephasings in the calculation. Taking the other factors into account, the total 
density matrix (ptot) of the photon pair is divided into three parts, 

1 



Ptot 



1 + 9 



(PPpol + (1 - P)pnoc + #Pnoise), 



(7) 



where 

Pnoc clXISGS from the FSS-induced distinguishability between the two decay paths, 
representing the non-overlapping part of the emission lines of excitons Xh and Xy, as 
shown in Fig.fTJ It describes those photon pairs which can be distinguished from the spectra 
domain, and have only a classical correlation without a phase relationship; hence, the non- 
diagonal elements of p noc are all zero, while the diagonal elements evolve in the same manner 
as ppoi- The third term p no ise, which describes the background noise, is set as an identity 
matrix. 
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FIG. 2: (a) Fidelity as a function of gate width w g {r g = 0) at S = 2.5 peV. Dashed horizontal 
line represents the limit for classical behavior, (b) and (c) present the fidelity as a function of r 9 
with gate width w g = 0.5 ns for S = 2.5 and S = 3.6 peV, respectively. 



In Ref. [13j, Stevenson et al. detected the photon pairs with delays r in the range r g < 
T < ( r g + w g ) by applying a single timing gate. In this case, the non-diagonal term p i4 
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(P41 = P14) has the form of an exponential decay 

Tg+Wg 

P14 oc J dr exp [— Tr/2 + iSr] 

p (iS-V/2)r g 

where T = 720 + 710 + T12 + 721 ■ Our results shown in Fig. [2] can explain their experimental 
behaviors as being a consequence of a maximum fidelity 0.73 of Bell state |$) obtained 
when gate width was 49 ps. Also, a marked oscillation in fidelity as a function of r g was 



observed, when the same parameters as that in Ref.[13j are chosen: 720 = 710 = 1-3 ns -1 , 
732 — 731 = 1-8 ns" 1 (that corresponds to the inverses of the exciton and biexciton lifetimes 
within a QD), and temperature T = 10 K, a value for r], the overlap part of the two exciton 
spectra, of 0.91 when 5* = 2.5 /zeV. To fit the experimental results, we chose g = 0.45. When 
we fix the gate width in Eq. ([8]), the fidelity as a function of r g has a sinusoidal envelop, as 
evident in Fig. 2 (b) for S = 2.5 fieV and (c) for S = 3.6 /ieV. The frequency of the oscillation 
increases as the FSS widens. Furthermore, we fix t 9 = 0, and change the gate width w g , in 
the long time scale, the integration of Eq. flH]) approaches zero. Consequently, the fidelity is 
below 0.5 (the maximum achievable fidelity for an unpolarized classical state). In contrast, 
when the gate width is shorter than 1 ns, the fidelity can increase above the classical limit, 
as shown in Fig.|2](a). 

However, fidelity is insufficient in describing the entanglement properties of the photon 



pair. To quantify this, a widely used form of entanglement is Wootters' concurrence 23], 
which is well defined. The concurrence C of a two-particle system is defined as 



C(p tot ) = max{0, V - V A 2 - V A 3 - a/ A 4 }, (9) 

where Aj are the eigenvalues of the matrix p to t(o"2 <8> cr 2 )pt t(o"2 <8> 02) arranged in decreasing 
order, and p tot is defined in Eq. ([7]) . 

The result is presented in Fig.[3] where the FSS energy is 2.5 /xeV. As temperature rises, 
concurrence decreases due to strong phonon induced scattering processes. In Fig.[3](a), 
entanglement declines rapidly for large gate widths, that can be ascribed to obtaining more 
which-path information from the FSS in the energy domain. This is because the Fourier 
transform of a truncated exponential decay results in a broad natural linewidth of the post- 
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selected photons. Long time transition processes between two exciton states also reduce 
entanglement, as shown in Fig.[3](b) where w g is fixed at 0.1ns. 

Under both situations, ESD is observed at high temperatures. It is worth noticing that 
this ESD is quite different from those in quantum optics experiments 6|, [7|, where one can 
alter the initial state to prevent the system from undergoing ESD. In a QD system, the ESD 
is mainly due to the intrinsic phonon assisted transition process between |1) and |2). Thus, 
we cannot change the initial state, but FSS can be tuned in experimnet, and we can show 
that the ESD phenomenon will disappear when the FSS is small and the temperature is low. 

Further, the influence of temperature and FSS are investigated, as shown in Fig.|H The 
photon pair generated by a QD with large FSS has lower degree of entanglement, that also 
falls off rapidly. However, for a QD with small FSS, such as S = 0.5 /ieV, the concurrence 
decays slightly with temperature, because a tiny phonon assisting emission rate results in a 
small transition rate in this small-FSS system. Comparing Fig.|H(a) and (c), for which the 
evolution time is the only differing parameter, entanglement at low temperature is found 
to be almost the same for every FSS, although the decreasing rates are apparently faster 
for longer evolution times (Fig.H](c)) and ESD temperatures are lower for larger FSSs. In 
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FIG. 3: (Color online). Concurrence as a function of temperature and (a) gate width w g (r 9 = 0), 
(b) delay time T g (w g = 0.1). The FSS is 2.5 fieV. 
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FIG. 4: (Color online). Concurrence as a function of temperature for different FSS energies. The 
four kinds of line shapes represent that S is 0.5, 2.5, 3.5, and 5 /ieV respectively, (a) w g = 0.1ns 
and T g = ns. (b) w g = 0.5 ns and r g = ns. (c) w g = 0.1 ns and r g = 0.5 ns. (d) w g = 0.5 ns and 
T g = 0.5 ns. 
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FIG. 5: (Color online). The relationship between the sudden death temperature and splitting S. 
Gate width w g = 0.1 ns. 

Fig.|H(b), the degree of entanglement is already small at low temperature due to its large 
gate width, and the reduction rate is faster than that in (a), a fact attributable to a system 
evolving within the time duration of the gate window. In Fig.|U(d), the influence of long 
evolution times and the gate widths are taken into account, so the concurrence initial values 
at low temperatures are small, and reduction rates are the fastest. 

Finally, we show that the observed ESD is independent of background light. In Fig. 5, we 
present the relationship of sudden death temperature as a function of splitting S for different 
g, for which we choose r g = 0.5 ns and w g = 0.1ns. Even in the ideal case, i.e. g = 0, we 
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still observe ESD. The ESD temperature decreases fast with respect to FSS, owning to the 
fact that for large FSS, not only is the transition rate between the two exciton states fast, 
but also the initial degree of entanglement is very small, due to the distinguished part p noc 
being dominant in the total density matrix p tot . This indicates that a large-FSS QD cannot 
generate good entangled photon pairs even at low temperatures. 

In conclusion, we have calculated the concurrence of photon pairs emitted from quantum 
dots and considered how temperature and fine structure splitting affect entanglement be- 
tween a photon pair. An enhancement of entanglement sudden death with temperature has 
been found. This obviously prevents quantum dots emitting entangled photon pairs at high 
temperatures. We also point out that this entanglement sudden death is a characteristic of 
quantum dot systems. 

This work was supported by National Fundamental Research Program, and National 
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